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HEAT KERNEL ESTIMATES FOR SUBORDINATE BROWNIAN 

MOTIONS 

ANTE MIMICA 


Abstract. In this article we study transition probabilities of a class of subor¬ 
dinate Brownian motions. Under mild assumptions on the Laplace exponent of 
the corresponding subordinator, sharp two sided estimates of the transition prob¬ 
ability are established. This approach, in particular, covers subordinators with 
Laplace exponents that vary regularly at infinity with index one, e.g. 

-1 or ^ P e (0.2) 


log(l -f A) 


log(l -f A/3/2) > 


that correspond to subordinate Brownian motions with scaling order that is not 
necessarily stricty between 0 and 2. These estimates are applied to estimate 
Green function (potential) of subordinate Brownian motion. We also prove the 
equivalence of the lower scaling condition of the Laplace exponent and the near 
diagonal upper estimate of the transition estimate. 


1. Introduction 

Recently, in probability theory and analysis there has been made progress in study of 
various properties of discontinuous Markov processes and their associated non-local 
generators. One of the notions that connects these two subjects is the heat kernel. 
In probability theory it can be understood as the transition density p{t, x, y) of a 
Markov process X, while in analysis it is the fundamental solution of the equation 
dtu = Au, where A is the inhnitesimal generator of X. Hence it is not surprising that 
one of the problems that drew much attention recently was to hnd sharp estimates 
of the heat kernel p(t, x, y) for various classes of discontinuous Markov processes X 
and non-local operators A. 

For pure jump symmetric processes with stable-like jumping kernels on or 
sharp heat kernel estimates were obtained in [BL02, CKOS, CKOS], while in [BGR14, 
Knl4, KS14, KS15, Sztll] they were obtained for some classes of Levy processes. A 
common property of all stochastic processes for which sharp two-sided estimates of 
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the heat kernel were obtained is that the scaling order was always strictly between 
0 and 2. This notion will be explained in detail later (see (1.2) and (1.3); for 
alternative approach to scaling Matuszewska indices may also be used, see [BGT87]). 
Our motivation was to obtain sharp heat kernel estimates of heat kernel when this 
property fails. Some far from optimal upper bounds in this case were obtained 
earlier in [Miml2] and [KS15]. 

Markov jump processes became important also in applications (e.g. in physics and 
hnance, see [CT]) and typical examples are a-stable processes, where a G (0, 2) and 
from this aspect it is also important to have good estimates of transition densities. 
Processes that were not covered by the theory known so far are conjugate geometric 
stable processes. These are Levy processes X = {Xt)t>o in such that, for some 
/?G(0,2], 


j- ICr 

The order of such processes is not strictly less than 2. Actually, concerning behavior 
of jumps and some other potential-theoretic notions (e.g. Green function or jumping 
kernel), these processes are between any rotationally invariant a-stable process {a G 
(0, 2)) and Brownian motion. It can be seen (see (1.6) and [Miml3, Miml2]) that the 
intensity of small jumps is higher than in the case of any stable process. Knowing 
sharp estimates of transition density of such processes might be useful in modeling 
various phenomena in nature by them. 

A very successful technique that was used to obtain upper bounds in heat kernel 
estimates was developed by Garlen, Kusuoka and Stroock in the paper [GKS87] and 
some of the already mentioned papers actually use this method. This method works 
well in stable-like cases (see [BGK09, BL02, GK08, GKK08, GKKll, KS12]), but 
it is not clear how to extend it to the cases when the scaling order is not strictly 
between 0 and 2. Hence, obtaining sharp heat kernel estimates for subordinate 
Brownian motion not satisfying this scaling order restriction could be a possible 
starting point for developing a generalized version of this method. 

Another motivation for this investigation was to try to generalize existing heat 
kernel estimates within the class of subordinate Brownian motions. As it will be 
seen from the main result, estimates will have a new form, but when the scaling 
order is restricted to be strictly between 0 and 2, these estimates reduce to the 
already known form. 

For a,b G M we denote a A b := min{a, 6} and a V b := max{a, 6}. Notation 
f{x) X g{x),x G / means that there exist constants ci,C 2 > 0 such that cif{x) < 
g{x) < C 2 g{x) for x G /. By Br{x) = {?/ G : |x — |/| < r} we denote open ball 
around x G with radius r > 0. We also use convention 0“^ = -|-oo. 
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A subordinator S = {St)t>o is an increasing Levy process, which is a stochastic 
process dehned on a probability space (f2, T, P) with stationary and independent 
increments with sample paths that are right continuous with left limits. It follows 
from the dehnition that S takes values in [0, cxo) and the Laplace transform of St is 
of the form 

A>0, 


where 0 is called the Laplace exponent of S. It has the following form (see [Ber96, 
IILl]) 

0(A) = b\+ [ (1 - e-^')p(df). 

J (0,oo) 

Here, 6 > 0 is called the drift of S and /i is a measure on (0, cxo) satisfying /(•qq^)(1 A 
t)fi{dt) < oo called the Levy measure of S . 


Let B = (i?o iPx)t>o,a;gR'i be the Brownian motion in {d > 1) independent of the 
subordinator S. We dehne the subordinate Brownian motion X = (W, Pa:)t>o,xGK'i 
by Xt = Bst, t >0. It is a Levy process (see [Sat99, Theorem 30.1]) such that 

Moreover, it has transition density p{t, x, y) = p{t, y — x) and it is of the form 

p{t,x) = f (47rs)“‘^^^e“^P(S't G ds). 

J (0,oo) 


Taking 0(A) = A“/^ with a G (0,2), we obtain rotationally invariant symmetric 
a-stable process X. Its inhnitesimal generator is the fractional Laplacian —(—A)“/^ 
dehned by 

-{-^T'^u{x) = [ {u{x + y)- u{x) - Xu{x) ■ |/lBi(o)(l/))^fe dy 
and the following heat kernel estimate holds (see [BG60]) 

X r^/“ A A (1.1) 

Note that in this case heat kernel estimate can be expressed just in terms of the 
Laplace exponent and its inverse function. Such type of estimate will continue 
to hold if the scaling of the process X is strictly between 0 and 2 (see [BGR14] 
and Gorollary 1.3), but if the scaling fails to satisfy this condition we will see that 
different form of heat kernel estimates appear. The approach in this paper is more 
general and it will essentially hold for subordinate Brownian motions with scaling 
order that is strictly greater than 0. 

We introduce the following scaling conditions for a function /: (0, cxo) —)■ (0, oo) : 
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(L) there exist 7 > 0 , Al > 0 

/(A) - 

(U) there exist 5 > 0, A[/ > 0 
/(At) 


/(A) 


< Cux^ 


and Cl > 0 such that 
for all A > Al and 
and Cu > 0 such that 
for all X > Xu and 


T > 1 , 


X > 1 . 


( 1 . 2 ) 


(1.3) 


If / is non-decreasing, then (L) and (U) are actually doubling conditions meaning 
that it is enough that they hold for some x > 1. 

Scaling conditions will be interesting for the Laplace exponent 0 of a subordinator 
S and for the function H : (0, 00) —)• [0, 00) dehned by H{X) := 0(A) — A0'(A). The 
function H appeared in the work of Jain and Pruitt [JP87], where, in particular, as¬ 
ymptotic properties of lower (upper) tail probabilities of subordinators were studied. 
We will need estimates of the tail probabilities of subordinators and in this sense 
tail estimates in Section 2 represent an upper tail counterpoint of the results from 
[JP87] (see Proposition 2.3 and Proposition 2.8) . However, our method involves 
different techniques. 

It is not hard to show that if H satishes (L) or (U), respectively, then the same 
holds for 0 in the case of a zero-drift subordinator (see Lemma 2.1). If 0 satishes 
(L), then processes with scaling order 0 cannot be considered, e.g. geometric stable 
processes. These are the processes obtained by subordinating Brownian motion by 
geometric stable subordinators, that is the subordinators with the Laplace exponent 
of the form 0(A) = log(l -|- A^/^) with 0 G (0,2]. The corresponding near diagonal 
estimate is inhnite for d> (3. Indeed, by [SSV06, (3.8)], 

poo /*|^|^ 

p(l,x) = / x) dt > Cl / df > C2|x|“'^+^, 

Jo Jo 

hence lim p(l,x) = J-cxd . Thus estimating heat kernel in this case makes sense 

1x1^0+ 

only away from the diagonal. 

Actually, the lower scaling property of 0 is equivalent to the near diagonal upper 
bound of the heat kernel and this is the hrst result of this paper. 

Theorem 1.1. Let S = (St)t>o be a subordinator and let X = (X0t>o be the corre¬ 
sponding subordinate Brownian motion in with the transition density p{t, x, y) = 
pit, y — x). Then there exist C > 0 and Al > 0 such that 

pit,x) < C(i)~^it~^Y^‘^ for all 0 < t < 0 (Al)“^ and x G (1.4) 

if and only if 0 satisfies (L)- 
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The main result of the paper is the following sharp heat kernel estimate. 

Theorem 1.2. Let S = {St)t>o be a suhordinator with zero drift and let X = {Xt)t>o 
be the corresponding subordinate Brownian motion in with the transition density 
p{t,x,y) =p{t,y-x). 

(i) If (f) satisfies (L), then there exists a constant k G (0,1) such that for all 

0 < t < and x G satisfying t(j){\x\~‘^) > 1, the following near 

diagonal estimate holds 

p{t,x) X (j)-\t-^Y/^. 

(ii) If H satisfies (L) and (U) with 6 < 2, then there exist constants G 

(0,1), C > 1 and ai, au > 0 such that for all 0 < t < and x G 

satisfying |x| < and t0(|x|“^) < 1, the following off-diagonal 

estimates hold 

p{fx) < C (t\x\-^H{\x\-^) V 
p{fx) > C-^ [t\x\-^H{\x\-^) V 


Note that, in the case \l = \jj = 0, the heat kernel estimates in Theorem 1.2 are 
global in space and time. Since f0(|x|“^) < 1 is equivalent to |xp0“^(t“^) > 1, we 
get the following off diagonal estimates too: 


p{fx) < C (^t|x|-‘'il(|x|-2) V 

p{fx) > C-^ (t\x\-^H{\x\-^) V , 


for some au > 0 and au = au- 

A novelty in this result is the appearance of the term of the form 

in the off-diagonal estimate. It could be explained as an intermediate term between 
on-diagonal estimate and the classical off-diagonal estimate t|x|“'^if(|x|“^) that in¬ 
volves tail estimate of the Levy measure of the subordinate Brownian motion X (see 
Figure 1) . 

A natural question that arises is in what situations the functions H and 0 are 
comparable. It turns out that this is equivalent to the property that 0 satisfies 
(U) with 5 < 1 (see Proposition 2.9). In other words, H can be replaced by 0 in 
Theorem 1.2 precisely when the scaling is strictly between 0 and 2 and in this case 
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t 



Figure 1 . Regions of heat kernel estimates. Dotted line corresponds 
to = 1, while ’fnll’ line corresponds to tH{\x\~‘^) = 

the estimate from Theorem 1.2 reduces to the estimate that is possible to deduce 
from results already known from previous works (see e.g. [BGR14, eq. (2)]). 

Corollary 1.3. Let S = {St)t>o be a subordinator with the Laplace exponent cf) satis¬ 
fying (L) and (U) with 5 < 1 and let X = {Xt)t>o be the corresponding subordinate 
Brownian motion in with the transition density p{t, x, y) = p{t, y — x) . Then 

there exist constants G (0,1) such that for all 0 < t < and x G 

satisfying |x| < l\p\y^‘^, 

p{t, x) X A . 

In some cases it is easier to work with the second derivative of 0 than H. 

Proposition 1.4. If (f is the Laplace exponent of a zero drift subordinator such that 
H satisfies (L), then there exists M > 1 such that 

X>MXl. 

A similar result has already appeared (see e.g. [JP87, Lemma 5.1]). 

Remark 1.5. Concerning comparability of certain function involving 0 there are 
more interesting results. It turns out that if H satishes (L) and (U) with 6 < 2, 
then the functions 0(A) and A0'(A) are comparable. Conversely, comparability of 
these two functions implies that 0 satishes (L) (see Proposition 2.10). 
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Let us return to the new features of our result and consider examples with order of 
scaling that is not strictly between 0 and 2. 


Example 1.6. Let 0(A) = 


0 -i(A) 


2 

A2-/3 

A log A 


iog(i+A/^/ 2 ) ^ where 0 G (0, 2). Here, 


0 < A < 2 
A > 2 


i/(A)x 


A1-/3/2 

A 

(log A) 2 


0 < A < 2 
A > 2. 


Hence, H satisfies (L) and (U) with A^ = A{/ = 0 and S < 2, but 0 and H are not 
comparable for large values of A (compare with Proposition 2.9). For example, if 
0 < t < 1/2 and 0 < |x| < 1/2 satisfy f0(|x|“^) < 1, then we get the off-diagonal 


estimate 


p{t,x) 


\x 


d+2 


(log R. 


Vf 


-d/2 


log- 


-d/2 


,nr 


where a = ai or a = au dependeing whether we consider lower or upper bound. For 
other t and x we can express estimates explicitly similarly by using behavior of H 
and 0“^. In particular, for t > 1/2 and |a;| > 1/2 we get stable-like estimates 


The hrst estimate gives an answer how a sharp upper bound in [Miml2] for a Levy 
processes should look like. In [KS15, Theorem 4] a similar upper estimate for this 
process has been obtained. 


Example 1.7. Let 0(A) = — 1 be the Laplace exponent of the conjugate 

gamma subordinator without killing. In this case 0(A) ~ -A and H{X) ~ ^ as 
A —)■ 0-I-. Hence, 


0-i(A) 


A 0 < A < 2 
A log A A > 2 


i7(A)x 



0 < A < 2 
A > 2. 


Therefore, 0 and H are not comparable at all, but H satishes (L) and (U) with 
Al = 0 and Xu = 2. In this case, we get the estimate only for small values of t and 
|a;|, which are basically the same as the corresponding estimates in Example 1.6. 


A crucial step in estimating heat kernel are upper tail estimates of subordinators 
(see Proposition 2.3, Proposition 2.5 and Lemma 2.6). Put together, if H satishes 
(U) with 6 < 2, we obtain the following estimate 

¥{St > r) X tH{r-^) (1.5) 

for t > 0 and 0 < r < MX^^ satisfying f0(r“^) < e“^(l — e), where M,e G 
(0,1) and constants in (1.5) depend on £. As already mentioned, estimate (1.5) 
could be understood as a uniform version of [JP87, Theorem 5.1] for upper tails of 
subordinators. 
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As an application, estimates of heat kernel obtained in Theorem 1.2 can be nsed to 
estimate Green fnnction (or potential) in the case of transient snbordinate Brownian 
motion X = {Xt)t>o- To be more precise, we say that X is transient if Po( lim \Xt\ = 

~ t^OO 

oo) = 1 . A necessary and snfficient condition for transience is that there exists a > 0 
so that jQ^^dy < oo (see [Sat99, Corollary 37.6]). If X is transient, then (see 
[Sat99, Theorem 35.4]) the following fnnction is well-dehned and hnite 

poo 

G{x,y) = G{x-y)= p{t,x-y)dt, x,yeR'^,x^y 

Jo 

and it is called the Green fnnction (or potential) of X . 


The following resnlt gives behavior of the Green fnnction when the scaling order 
is strictly greater than 0. Snch estimates are usnally proved by using different 
techniques (e.g. using potential measure of the subordinator), but since we have 
obtained sharp heat kernel estimates it is almost a direct consequence. 


Corollary 1.8. Let S = {St)t>o be a subordinator with zero drift such that its Laplace 
exponent (f satisfies (L) with Al = 0 and let X = (Xt)t>o be the corresponding 
subordinate Brownian motion in Rfi. Assume that X is transient and that in the 
case d < 2 the Laplace exponent satisfies additionally (U) with 6 < d/2 and 
Xu = 0. Then 


G{x) 


1 

x\‘^(j){\x\~‘^) ’ 


X . 


Green function of the conjugate gamma process from Example 1.7, for d > 3, can 
be estimated as 


G{x) 


\x 


\x 


2-d 


2-d 


log 


hi 


X < 


X 




( 1 . 6 ) 


This shows how close this process is to the Brownian motion, since the Green func¬ 
tion of the Brownian motion in (d > 3) is the Newtonian potential G^^)(x) = 


c X 


2-d 


Estimate in Corollary 1.8 for subordinate Brownian motions is known in many cases 
when the scaling order is strictly between 0 and 2 (see [SV09, KSV12]). For some 
cases it has been calculated when the scaling order is allowed to be near 2 (see 
[Miml3]) and recently it has been proved that the condition (L) for is actually 
equivalent to the lower bound of the Green function in Corollary 1.8 for |x| small 
(see [Grzl4]). The form of the Green function behavior in the case of zero order 
scaling is different (for details, see [SSV06, KM12, Miml2]) . 

The structure of the paper is as follows. In Section 2 tail estimates of subordinators 
are obtained. Further, some properties of functions and H and their relationships 




HEAT KERNEL ESTIMATES FOR SBM 


9 


are proved. Some of the results in this section could be also considered to be of 
independent interest. Section 3 is devoted to heat kernel estimates. It consists of 
three subsections, hrst two for upper and lower estimates and in the last one we 
prove Green function estimates and the main results. 

2. SUBORDINATORS AND THEIR TAIL ESTIMATES 

Let S = {St)t>o be a subordinator. Recall that the Laplace transform of St is given 
by 

E[e-^^‘] = 0(A) =6A+/' (1 - e-^^)/i(d|/), A > 0, 

J (0,oo) 

where 0 is called the Laplace exponent, 6 > 0 is the drift and /x is the Levy measure 
of S meaning that /x is a measure on (0, oo) satisfying ^^(1 A y)fi{dy) < oo . 

The Laplace exponent 0 belongs to the class of Bernstein functions BT = {/ G 
G°°(0,oo): / > 0, (—> 0, n G N} . It is known that every Bernstein 
function / has a unique representation 

/(A)=a + 6A+/ il-e-^y)y{dy), (2.1) 

J (0,oo) 

where a, 6 > 0 and /x is a Levy measure (cf. [SSV12, Chapter 3]). 

The following formula that follows from Fubini theorem will be useful 

= — + 6 + [ e~^^y{t,oo) dt, A>0. (2.2) 

i(0,oo) 

Further, the following inequality will be often used: 

A/'(A) = Xb+ [ Xte-^^y{dt) 

J (0,oo) 

<Xb+ [ (1-e“^‘)/x(dt) </(A), A>0. (2.3) 

J (0,oo) 

A direct consequence of inequality (2.3) is non-negativity of the function H. 

The hrst lemma gives some general properties of Bernstein function needed in the 
sequel. Although most of results are known, we present some alternative proofs. 

Lemma 2.1. Let 0 G BJ^. 

(a) For any A > 0 and x > 1, 

4>{Xx) < x(j){X) and H{Xx) < x^H{X). 
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(b) Assume that the drift off) in the representation (2.1) is b = 0. If H satisfies 
(L) (resp. (U) then f) satisfies (L)(resp. (U) 


Proof, (a) Without loss of generality we may assume that a = 6 = 0 in the represt- 
ntation (2.1) of 0. The result concerning 0 is already known and follows from 
concavity of 0. Here we present new proofs of both inequalities that follow from the 
following elementary inequalities 

1 - < x{l - e ~^), 

1 - e-^^ - \xe-^^ < x2(l - e-^ - Ae”^), x > 1, A > 0 . 


Indeed, 


0(Ax) = / (1 — e ^)fi,{dt)<x / (1 — e )iJ.{dt) = x(j){X) 


’ (OjOo) 


( 0 ,oo) 


and 


H{Xx) = (1 - - Xxte-^^^)fx{dt) 

J (0,oo) 


< X 


(1 - - Xte-^^)fi{dt) = x^H{X) 


.-At 


' (0,oo) 


(b) Assume now that H satishes (U) and that the drift is 6 = 0. Then lim = 
6 = 0 (see [Ber96, Proposition 1.2 (ii)]); hence 


A^oo 


X 

By change of variable, 

poo 

0(Ax) = Xx 

< CuXx^ 

Similar argument applies to (L). 


H{s) 


ds. 




H{sx) 


ds 


Xx s 

H{s) 


'X 


ds = Cux^f){X). 


□ 


Note that Lemma 2.1 (a) suggests that 0 and H could have different scaling proper¬ 
ties and this can be seen in Example 1.7. To be more precise, for 0(A) = — 1 

the bounds from Lemma 2.1 (a) are attained: 


0(Ax) 

hm 

A—>-0+ 


= X 


and 


lim 

A^0+ H{X) 


The following tail estimate of the subordinator reveals a probabilistic connection 
between the functions H and 0 and it will play an important role in obtaining 










HEAT KERNEL ESTIMATES FOR SBM 


11 


upper off-diagonal estimates of the heat kernel. Before this we need an auxiliary 
lemma. 


Lemma 2.2. Let /i he the Levy measure of a suhordinator with the Laplace exponent 
(f and the Levy measure p. Then, for any r > 0, 


(i) p{r,oo)+r / yp{dy) = 2r / yp{y,oo)dy 


(0,r] 


(ii) / yn{y,oo) dy < er^H{r ) 


(in) / yp{dy) < ercfir 




'{0,r] 


Proof, (i) By Fubini theorem we get 


' (0,r] 


y‘^ti{dy) = 


'{0,r] Jo 


2zdzpi{dy)= / 2zpi{{z,r]) dz 


= / 2z{jj,{z, oo) — jj,{r, oo)) dz = 2 / zii{z,oo) dz — r'^p,{r,oo). 


(ii) Taking derivative in 

it follows that 

H{\) 0(A) - A0'(A) 


0(A) 


A2 


A2 


= b+ e ^^niy, oo) dy 


e ^yyn{y,oo)dy> e ^ yp,{y,oo)dy. (2.4) 


The desired inequality is obtained by choosing A = r ^. 

(iii) Since 1 — e~^ > xe~^ for t > 0, 

0(A) > [ \ye~^yp,{dy) > e~^\ [ yp{dy) 
J(0,oo) >''( 0 ,A“h 

and now it is enough to take A = r 


-1 


□ 


Proposition 2.3. Let S = {St)t>o be a suhordinator with the Laplace exponent 0. 
Then there exists a constant C 5 > 0 such that 

lP(5'i > r(l -I- et0(r“^))) < CstH{r~^) 

for all t,r > 0 . In particular, if 0 < e < 1, then 

^{St >r)< CstH{e-^r-^) < Cse-HH{r-^) 

for all t,r > 0 satisfying t(l){r~^) < e“^(l — e). 
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Proof. Let 

St:=St-tl yfi{dy), t>0,r>0. 

J (0,r] 

Then 


^ ^ / yfi{dy) 

J{0,r] 


= exp <-t (l - + iOyl{\y\<r}) y{dy) 

J(0,oo) 


which shows that {5't}t>o is a Levy process in R with the Levy exponent fjiO) = 
/( 0 oo)(^ — + i9yl{\y\<r})p{dy). Hence it follows that for any / G Co(®) 

{/’e ^^(R): e C'o(R)}, where C'o(R) = {/ G C'(R) : lim /(x) = 0}, the 

inhnitesimal generator of S is given by (see [Sat99, Theorem 31.5]) 


^f{x)= / {f{x + y)-f{x)-f'{x)yl{\y\<r}) y{dy). 

J (0,co) 

Let g G C^(R) with bounded hrst and second derivatives such that 

1 3 

0 < g{x) < 1 for all X G R, g{x) = 0 for x < - and g{x) = 1 for x > - 

and dehne, for n > 2, 

{ g{x) X < 1 

1 1 < X < n 

g{n + 1 — x) X > n . 

Then /„ G C'q(R) and 

ll/nlloo = 1, ll/nlloo < Us'lloo aud ||/"|loo < ll/IU for all neN. 

Let n G N and set fnAv) ■= fn^) ■ Then 

fn,rA + y)- fn,rA) “ -^fnAAV^llvl^r} = fn (^) “ fn (f) “ A (f) 

— 2ll/nll°° (r) fo ll/^^l|ool{|^/|>r} 
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and Lemma 2.2 yields 



CsH{r-^) 


with Cs '■= 2e(||5'"||oo V 1). By Dynkin formula (see [RY05, Proposition VII.1.2 and 
Proposition VII. 1.6]) and the last display it follows that 



Jo 

Since t lim fn,r{x) = gr{x) := g{^) we may use monotone convergence theorem to 


conclude that 


= lim nUASt)] < CstH{r-^) . 


Hence, Lemma 2.2 (hi) and the last display dually yield 


P(^t > r(l + et0(r-i)) < F{St > r) < E[^,(^i)] 



Jo 

Let 9 > 0 and t,r > 0 such that t0(r“^) < 9 , then 

ASt > r(l + e9)) < ¥{St > r(l + etAr~^)) < CstH{r-^). (2.5) 


Now, if tAr ^) < by Lemma 2.1 (a), 

f0(r“^(l + e9)) < f0(r“^)(l + e9) < 9, 

hence, by (2.5) and Lemma 2.1 (a), 

P(5t >r)= F{St > r(l + e9)-^l + e9)) < Cstff(r-\l + e9)) 

< (1 + e9fCstff(r -^). 

Now take e = -As ■ D 

l-hetj 

The following estimate will be needed to get the lower bound estimates of heat 
kernel. 

Proposition 2.4. Let S = {St)t>o be a suhordinator with Laplace exponent 0. For 
any a, /3 > 0 such that ajd > 1 there exist constants p G (0,1) and r > 0 such that 



for all t > 0 . 
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Proof. Let ajS > 1. Since y ^ 1 — e ^ is strictly increasing and y ^ e ^ is strictly 
decreasing, Markov inequality implies 


P 


o‘4> H/3 L 


4-TT < ) > 1 -P ) -P < 


0 Mp* L 


acp 


= 1 - P (^1 - e"'^ > 1 - - P (e-'f’ 

_ p-0“dpi“L‘S'tl F L-0“L/3*“L>s't 


E 


> 1 - 

= 1 - 
= 1 - 


1 - e-^ e-L“ 

X _ g-i0(0“dpi“L) g-i0(0“d/3t“D 


1 — e 


-1 


3 -l/a 


3-P 


_ gl/o-/5 


1 — e“^ 

Since 1 — > 0 we can choose p G (0,1) small enough so that r := 1 — — 

gl/a-/3 > Q _ □ 


To obtain the lower bound the following result will play an important role. 

Proposition 2.5. If S is a subordinator with the Levy measure p and the drift term 
b > 0, then for all t,r > 0 we have 

> r) > 1 - . 

Proof. Let r > 0 and let T = (Tt)t>o and V = (fd)t>o be two independent Levy 
processes with Levy measures 

p ■ 1 ( 0 ,r] and p ■ l(r,oo) respectively , 

and T having drift b. Then St = Tt+Vt for t > 0 and T and V are both subordinators. 
Note that P is a compound Poisson process (since its Levy measure is hnite). Then 
we have 

F{St >r)> F{Vt > r) = 1 - 

since the event {Id > r} occurs precisely when the the compound Poisson process 
V jumps for the hrst time before time t. 


□ 

Since the tail of the Levy measure appears in the estimate in the previous result, it 
would be useful to know its behavior. 

Lemma 2.6. Let p be the Levy measure of a subordinator with the Laplaee exponent 
(j). Then 

p{r\ oo) < 2eH{r~^) , r > 0 . 
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If H satisfies (U) with 6 <2, then there exists M G (0,1) such that 
/i(r, oo) > 0 < r < 

Remark 2.7. Note that for A;/ = 0 the estimates hold for all r > 0. If H satishes 
(U) with 5 < 2, we may write 

/i(r, cxd) X H{r~^), 0 < r < MX^^. 


Proof. Taking derivative in (see (2.2)) 

0(A) 


A 


we get 


H{X) 

A2 


= b+ e oo) dy 


e ^yyy{y,oo)dy, A > 0 


( 2 . 6 ) 


and changing variable it follows that 

poo 

H{r~^) = / e~^yy{ry, oo) dy, r>0. 

Jo 

The npper estimate follows immediately, since it follows from (2.7) that 

pi 

-1 


(2.7) 


hence 


H{r )> e ^yyiry, oo) dy > ^y{r, oo), 


jj.{r, oo) < 2eH {r ), r>0. 


( 2 . 8 ) 


To get the lower bonnd we start with (2.7) and nse the npper estimate (2.8) and 
(U) to get, for any M G (0,1) and 0 < r < A^^, the following 

pM poo 

H{r-^)= e~yyy{ry,oo)dy+ e~yy[i{ry,oo) dy 
Jo Jm 

pM poo 

<2e e~^yH{r~^y~0 dy + y{Mr,oo) / e~^ydy 


' M 


< 2eCu{2 - 6)-^M^-^H{r-fi + (1 + M)e-^y{Mr, oo). 

Choosing M G (0,1) so that 2eC';7(2 — 6)~^M'^~^ < A and nsing Lemma 2.1 (a), for 


0 < r < A[j , we get 


y{Mr, oo) > 


2{1+M) 


H{r-fi > 


2(1+M) 


H{{Mr 




It is enongh to note that — "T 


□ 
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Proposition 2.8. Let S = {St)t>o be a zero-drift subordinator such that H satisfies 
(L) and (U) with 6 < 2. There exist constants M G (0,1), L > 1 and C 5 > 0 such 
that for all t > 0 and 0 < r < A MXfj^ satisfying t(j){r~^) < 1, 

P(r < St < Lr) > cstH{r~^) . 


Proof. By Lemma 2.1, Proposition 2.3, Proposition 2.5 and Lemma 2.6 it follows 
that there exists M G (0,1) such that, for r G (0, A MX'jf^), 

P(r <St< Lr) > ¥{St > r) - F{St > Lr) 

> 1 - - 4CstH{L-^r-^) 

> tfi{r, - 4CsCf^L-^tH{r-^) 

J\/f2 

> —- 4CsCf^L-nH{r-^) 

i\/T2 

> —- 4CsCf^L-pH{r-^) 

> tH{r-^) f —e-'" - 4CsCf^L-A 


In the first inequality we have used that 

< CE'L-'t4,{r-^) < Cin-i < L 

for L > 1 chosen large enough, since (L) also holds for 0 by Lemma 2.1. Again, 
choosing L > 1 large enough so that ^e~^^ — 4CsCf^L~'^ > 0 we obtain the desired 
bound. □ 

Now it is possible to prove the result that establishes comparability of Lf and 

xH-f>"(x)). 


Proof of Proposition l.f. Here we use elementary inequalities 


\x^e ^<l-e ^-xe x > 0. 

Then 

iA2(-f'(A))= [ |AVe-^V(^^^) 

J (OjOo) 

< [ (1 - - Xte-^^)pi{dt) = H{X), A > 0 . 

J (0,oo) 


HEAT KERNEL ESTIMATES FOR SBM 


17 


For M G (0,1) we obtain 

H{X) < f (1 — + n{M~^X~^ , oo) 

< [ XH^e-^^fi{dt) + 2eH{MX) 

< Xe^~"X\-(j)''{X)) + 2eCZ^M"'H{X), 

where in the second inequality we have used Lemma 2.6 and in the last (L). Choosing 
M G (0,1) so that 2eCj^^M'^ = | we get 

H{X) < e^''A2(-0"(A)), A > M’^Al . 


□ 

This section ends with results that explore comparability of the functions H, 0 and 
A0'(A). 

Proposition 2.9. The upper scaling condition (U) for (j) with 5 < 1 holds if and 
only if 

/f(A)x0(A), A>At;. 


Proof. Assume first that (U) holds for 0 with h < 1. Changing variables in (2.6) 
we obtain 

/•OO 

4>{X) = / e~'^n{X~^y, oo) dy, A > 0 . 


Since H{X) < 0(A), by Lemma 2.6, for any M G (0,1) we have 

rM 


noo 


0(A) = / e V(A ^y,oo)dy+ / e ^/i(A ^y,oo)dy 


rM 


' M 

, -1 


<2e 4>{Xy )dy + iJ,{MX ,oo) / e ^ dy 

Jo J M 

i-M 

<2eCu<f(X) / y-^dy + 2eH{M-^X)e-^ 

Jo 

< 2eCu{l - 6)-^M^-^(j){X) + 2eM-^H{X), 


where in the last inequality we have used Lemma 2.1 (a). Since 1 — h > 0, we can 
choose M > 0 small enough so that 2eCu{l — = | fo g^t that 


^0(A) < H{X) < 0(A) for all X > Xu 
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Assume now that there exists c G (0,1) such that 0(A) > H{X) > for all 

\> Xjj. This implies 


0'(t) ^ 1 - C 


t > X 


u ■ 


t 

Integrating this inequality over [A, Ax] for X > Xu and x > 1 we obtain 

0(Ax) 


log 


0(A) 


< (1 — c) logx 


which is (U) with 5 = 1 — c<l. □ 

The last result of the section explores comparability of functions and lower scaling 
conditions. 

Proposition 2.10. The lower and upper scaling conditions (L) and (U) for H with 
6 < 2 imply 

A0'(A)x0(A), A>AiVM-^At7 , (2.9) 

for some constant M G (0,1) . Conversely, 

A0'(A) X 0(A), A > Ai 
implies the lower scaling condition (L) for 0 . 


Proof. The upper bound follows from (2.3). Assume that (L) and (U) hold for H 
with 5 < 2. Then, for L > 1, Lemma 2.6 yields 
/•A-i 

A0'(A) > / Xte-^^fr{dt) > {fr{L-^X-\ oo) - /i(A“\ oo)) 

> L-^e-^ {c,H{XL) - C 2 H{X)) = c^P-^e-^HiX) 

> (A) (ciC^^C'lL'’' - l) , A > M-^Xu 

for some M G (0,1). Choosing T > 1 large enough so that — 1 > 0 we 

get A0'(A) > c^H{X) for all A > A^ V M~^Xu implying 

A0'(A) > ^^0(A) for all A>AlVM-^Ac/. 

1 + C3 


Assume that there exists c G (0,1) so that A0'(A) > c0(A) for all A > A^. Similarly 
as in the proof of Proposition 2.9 we obtain that for all x > 1 and X > Xu the 
following inequality holds 


log 


0(Ax) 

w 


> clogx 


and this is (L) with 7 = c. 


□ 
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3. Heat kernel estimates for SBM 

In this section we obtain estimates of transition density of the subordinate Brownian 
motion. 

Let S = {St)t>o be a subordinator defined on a probability space P) with 

the Laplace exponent 0 and let B = Pa.)^>o^ 3 ;eRd be the Brownian motion in 

> 1) independent of S. The subordinate Brownian motion is a stochastic 
process X = (X*, Pa;)t>o,rrGRrf defined by X^ = f > 0. 

Recall that it is a Levy process with the characteristic (Levy) exponent (/'(I'CP) and 
it has transition density given by 

p{t,x,y) = p{t,y - x) = [ e ds) (3.1) 

J (0,oo) 

for x, y G and t > 0. 

The following observation is important in obtaining estimates of transition density. 
Since B and S are independent, we may rewrite Po(|Xt| > r), for f,r > 0, as 

Pod^il >r)= j Po(|B,d > r2)P(,St G ds) 

J (0,oo) 

= [ Fo{{2s)-^\B,\^>{2s)-^r^)F{Steds) 

J (OjOo) 


Using the fact that chi-square distribution with parameter d (as a sum of 

squares of d independent standard normal random variables) and since it is inde¬ 
pendent of S, we conclude that 


where 


and T(t) 


Pod^tl >r)= Fo{Y> (2s)-V2) F{St G ds) 


Po(>^ >t) = 


' (0,oo) 

f 

qo,oo) 

1 


2'^/2r(d/2) 


nSt > (2|/)-V^)Po(X G dy ), 


gdl2-l^-sl2 


J^y* ydy, f > 0 is the gamma function. 


(3.2) 

(3.3) 


3.1. Upper bounds. We start with a lemma that shows that the inverse function 
of (j) also satisfies certain scaling properties. 

Lemma 3.1. Let (p he a Bernstein funetion. 
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(i) Let I := lim 0(A) and u := lim 0(A). Then for any I < \ < u and x > 1 

A—>0+ A—>-00 

such that Xx < u we have > x. 

(ii) If (L) holds for 0, then 


0 ^(Ax) 

TW 




for all A > 0(Al), x > Cl ■ 


Proof, (i) Since 0(Ax) < 0(A)x for any Bernstein function 0 (see Lemma 2.1), A > 0 
and X > 1, we can rewrite this (using the fact that 0 is strictly increasing) as 

0~H0(A)x) ^ 

0-H0(A)) 

Setting T] = 0(A) we obtain > x for all I < r] < u such that rjx < u. 

(ii) Let X > Xl and x > 1. Since 0“^ is increasing, one can rewrite (L) as 


0-HClx^0(A)) ^ 
0-H0(A)) 


Taking rj = 0(A) and y = Clx'^ it follows that 


0 Hw) 

0 - 1 ( 77 ) 


<Cl 


“I/o',,1/7 


y 


r] > 0(Al), yPCi. 


□ 

We proceed with the proof of the near diagonal upper bound for the transition 
density. This type of estimate can be found in [KS12] for more general Levy processes 
via Nash inequality. Our estimate is obtained by using scaling properties of the Levy 
exponent. 

Proposition 3.2. Assume that (L) holds for 0. There exists a constant Oi > 0 
such that 

p{t,x) < Cl 

for alio < t < (P{Xl)~^ and x G . 

Proof. Since e^^'^p{t, x) dx, the Fourier inversion formula yields 

p(t,x) = (27r)-'^ / • (3-4) 

jRrf JRd- 
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Switching to polar coordinates we get 

pCXD poo 

Pit, x) < Cl / dr = Cl / 

Jo Jo 

/•oo 


e ^ ds dr 




= Cl 


drds = f 


e-^ [r\j)]^ ds 


Now we split the integral in the last display and use monotonicity of 0 ^ and Lemma 
3.1 again to obtain 

/■l j roo . 


pit, a:) < t I / <?”['#> ‘(f)] '‘ds+ e * p ‘(f)] ’ ds 


< t [T‘(f)] 

< t [r‘(f)] 


e ds + 


^-hf) 


ds 


d r°° 


e^ds + Cj^'' e ds 


It is enough to take Ci ■= ^ ) Jq ^ ^ ds + Cj^ e ds 


□ 


The next result is the off-diagonal upper bound. We start with an estimate of the 
tail of the chi-square distribution given by (3.3) . 

Lemma 3.3. For any to > 0 there exists a constant c = c(to) > 1 such that 
^_Yi2-i^-ti2 < >t) < t>to. 


Proof. Using integration by parts, 

/ OO poo 

= -Vid-2) //2-2e-/2ds. (3.5) 

Assume hrst that d > 2. Since ds < t~^f{t ), we get /(t)(l — t~^d) < 

2 ^d/ 2 -ig-t/ 2 ^ yielding 

/(t) < t>2d. 

Therefore, in this case, from (3.5) we conclude that, for some constant ci = ci(to) > 

4, 

< fit) < t > to • 

Let d = 1. By (3.5) we get /(t) < . On the other hand, 

/ 2t 

S-^/\-^/^ds > 2^/2t-l/2(g-U2 _ g-i) > 21/2(1 _ g-io/2)^-l/2g-t/2 

for t > to . □ 
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Proposition 3.4. Assume that (L) holds for (p. There exist constants C 2 > 0 and 
Of/ > 0 such that 


p{t, x) < C 2 




for alio <t < (2e0(Ai)) ^ and a; G satisfying t(p{\x\ ^) < 1 . 


Proof. Let r > 0, 0 < t < (2e0(Ai)) ^ be such that t(j){r < 1. We first note that 
r‘^(j)-\t-^) > 1 and 2e << (2e)^/^ , (3.6) 

where the second inequality follows from Lemma 3.1 . Let B = be (d + 2)- 

dimensional Brownian motion and denote by X = (Xf)t>o corresponding subordi¬ 
nate Brownian motion. Let X be a random variable with the chi-square distribution 
with parameter d -|- 2. Using (3.2) and Proposition 2.3 it follows that 


P(|Xt| >r) < 




lt(f>(2yr ^)>7^ 


ACstH{2yr-^)Fo{Y G dy) 
Po(? G dy) 


■ 2e 

< 4CstH{r-^) [ (1 + {2yf)Fo{Y G dy) 

J (0,oo) 

+Po (y >yr\ht-^)). 


since 


< 1 V (2i/)2 by Lemma 2.1 . 


(3.7) 


Using Lemma 3.3 we obtain 

P(|^t| > r) < Cl (^tid(r-2) ^ 

since (2e)~^'^^ by (3.6) . 

On the other hand, from (3.1), 

Pd^il >r) = C2 f f P(^t G ds) dx 

JJ (0,oo) 

= C3 [ f dyF{St e ds). 

J (0,oo) J r 


(3.8) 
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Integration by parts yields 


' 2 2 
yd+i^-^ _ 2sr‘^e~i^ + 2ds 


d-l 

y e 4 s 


dy 


r! 

> 2sr e 4 s 


Hence, 

IP(|-^t| > r) > 2 c 2 ,r^ / G ds) = 2 c 3 |a;|'^p(t, t) , 

J (0,oo) 

for a; G such that |a;| = r by (3.1) . Using the last display (3.8) and (3.6) we get, 
for some Of/ > 0 , 

|a;|‘^p(t,a:) < C 4 (tH{\x\-^) + |a:|'^ 0 -^(r^)'^/ 2 g-aa|xp<A-i(t-i)j 
0 < t < {2e(f){\L))~^ and a; G satisfying t0(|a;|“^) < 1. 


□ 

3.2. Lower bounds. In this subsection the lower bounds of the transition density 
will be proved. 

Proposition 3.5. Assume that (L) holds for (p. There exist constants 6*4 > 0 and 
n G (0,1) such that 

p{t,x) > C' 4 (/>"^(f^) 5 e- 2 N''^“'(*'') 

for all 0 < t < k(1){\l)~^ and a: G . In particular, if additionally t0(|x|“^) > 1 
holds, we have 

p{t,x) > C 50 "^(t"^)^ , 

where = C 4 e“^ . 

Proof. First we note that (L) implies lim ^(A) = cxo. Let 0 < t < 0(Ai)“^ and 

A—>-oo 

X G . Using (3.1) we get 
p{t,x) = ci f G ds) 

J (0,oo) 

> Cl / G ds) 

> ci0-^(pr^)'^/2g-ip|2<^-irLp (2-V^(f^)“^ <St< 0-^(pr^)-^), (3.9) 

where a = 2, fi = 1 and p G (0,1) are the constants from Proposition 2.4. By 
Lemma 3.1, for 0 < t < p0(Ai)“^, 

rH.pt-') = r'(r>)£^ > ci'v'-'rHn. 
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Using the last display and Proposition 2.4 in (3.9) we get 

If t(f){\x\~‘^) > 1, then \x\^< 1. Hence, from the last display we conclude 

p{t,x) > . 

□ 

Proposition 3.6. Assume that (L) and (U) hold for H with 6 < 2 and that the 
drift of the subordinator is zero. There exist constants Cq > 0, L > 1 and M G (0,1) 

such that for all 0 < t < and x G satisfying |x| < 

and tcj){\x\~‘^) < 1 we have 

p{t,x) > CQt\x\~'^H{\x\~'^). 

Proof. From Proposition 2.8 it follows that there exist M G (0,1), L > 1 and > 0 
such that for |xp < L A MXj/ it holds that 

P(|x|^ < St < F|x|^) > cstH{\x\~‘^). 

Hence, by (3.1) we obtain 

p{t,x) > Cl f s“'^'^^e“'^P(S'i G ds) 

J[\x\^,L\x\‘^] 

> ciL-^/^\x\-'^e-^/^¥{\x\^ <St< L\x\^) 

> C2t\x\~^H{\x\~‘^) . 


□ 


3.3. Proofs of main results. In this section we prove our main results. We start 
with the proof of the equivalence of the lower scaling condition (L) for 0 and the 
near diagonal upper bound of the heat kernel. 

Proof of Theorem 1.1. If (L) holds for 0, then we use Propositon 3.2. Assume that 
(1.4) holds. Then 

[ {A7rs)-’^/^¥{St G ds) = p{t, 0) < 

J (0,oo) 

for 0 < t < 0(A£,)“^. In particular, for a > 1 we have 

C'i 0 -i(r ^)'^/2 > f {47rs)-^/^¥{St G ds) 

J (O,2a0“^ 

> (87ra)“'^'^^0“^(af“^)'^/^P(5i < 2a0“^(af“^)“^). 
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Since 


F{St < 2a(p-\at-^)-^) = 1 - P (1 - > 1 - 

_ g-i0(0“^(ai“^)) g-a _ g-2a 


> 1 


1 — e 


—2a 


1 — e 


—2a 


=: Cl > 0 


we obtain 


< 87ra(C'/ci)^/'^ =: 6, 0 < t < 0 (Al) ^ 

0 1) 


Note that 6 > 1. Taking A = 0 in the last display we get 

0(A6) 


0(A) 


> a, \ > \l ■ 


(3.10) 


Let X > 1 and choose n G N so that 6 ” ^ < x < If’. By iterating (3.10), for any 
A > Al we get 


0(Ax) ^ 0(A6"' ^) ^ . 

- ~wr - “ 



Hence, (L) holds with 7 = 


□ 


Proof of the main result is a consequence of propositions established in the previous 
subsections. 


Proof of Theorem 1.2. Follows directly from Proposition 3.2, Proposition 3.5, Propo¬ 
sition 3.4 and Proposition 3.6 since (L) also holds for 0 by Lemma 2.1. □ 


Proof of Corollary 1.3. Note hrst that, by Proposition 2.9, 0 and H are comparable, 
so (L) (resp. (U)) are basically the same for these two functions. 


Assume that t0(|x| ^) < 1. Then |xp0 ^{t ^) > 1 and so there exists a constant 
Cl > 1 so that 

Hence, by Lemma 3.1 (i). 


0 {t 


-Cd/2^-au\x\ 


2^-lU-l 


(i L < 


S Cl X 


-d 




< Ci|x| '^t0(|x| ^), (3.11) 


implying the upper bound 
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On the other hand, by Proposition 2.9, 

for \x\~‘^ > Xu . The estimate follows now from Theorem 1.2 . □ 

It is left to prove the Green fnnction estimates using the heat kernel estimates. 


Proof of Corollary 1.8. We split the integral 

G{x) = / p{t,x)dt+ / p{t,x)dt 

Jo J 

and use already established upper and lower bounds. By Proposition 3.5 and Lemma 
3.1, 

POO ^2(^(P|-2)-i 

/ p(t,x)dt>ci / (j)~^dt 


2 J (j){\x\~'^) \x\^(j){\x\~‘^) 

On the other hand, using the near diagonal upper bound from Proposition 3.2 we 
get 




p{t,x)dt<c2 / 0 ^Y^‘^dt = c2 I Jy 


= C2 


G(h|-2)-i 

P\x\-^ 

/ y 


dl2 

0 V 0(l/) 


JO J^iyy Jo (l>iy) 

since XcfJX) < 0(A). If d > 3 we take h = 1 by Lemma 2.1 and for d = 1, 2 we use 
(U) and the assumption 26 < d. In both cases we have | — 1 — d > —1. Hence 


'hi ^ . 412-1 


^ dy< 


’hi’ 


,d/2-l 


\X 


2 

-2N <5 


dy = 


C4 


Jo 4>{y) 0(kl"^) Jo \ y J |x|'^0(|a:|-2) 

It is left to estimate the integral of the off-diagonal upper bound from Proposition 
3.4 integrated over the appropriate time-interval: 

/ p{t,x)dt<c^ U|x|-‘^id(|a;|-2)+ ^ )\ dt 

< C5(|x|“'^id(|a:|“^)-|-C4|x|“'^0(|a;|“^)) / tdt 


1 , , ^,'H{\x\-^) 1 

= |c5 X -'^ ( ,T + 


< 


C5 


^ ' 0(|x|“2)2 0(|x|“2) J |x|'^0(|x|“2) ’ 

where in the second inequality we have used (3.11) and in the last H{X) < 0(A). 
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